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Abstract. We report on a study of the pentaquark ©^(1540), using a variety of different interpolating 
fields. We use Chirally Improved fermions in combination with Jacobi smeared quark sources to improve 
the signal and get reliable results even for small quark masses. The results of our quenched calculations, 
which have been done on a 12^ x 24 lattice with a lattice spacing of a = 0.148 fm, do not provide any 
evidence for the existence of a O'^ with positive parity. We do observe, however, a signal compatible with 
nucleon-kaon scattering state. For the negative parity the results are inconclusive, due to the potential 
mixture with nucleon-kaon and A^*-kaon scattering states. 

PACS. 12.38.Gc Lattice QCD calculations 



1 Introduction 

The possible discovery of the 6»+(1540) by the LEPS Col- 
laboration at SPring-8 Tl has initiated great interest in 
exotic baryons. Since then, there has been a large num- 
ber of experiments that have confirmed this result P] 
but also about the same number that could not confirm 
it 23]-|2ni- Presently the experimental situation is highly 
unsettled |12P^- 

Also lattice QCD has not been able to give a conclusive 
answer so far to the question whether there might exist a 
bound 0+ state or not 

To confirm or disprove the existence of the 0+ by lattice 
calculations is a huge challenge. If a clear conclusion of the 
could be reached and would eventually be experimentally 
confirmed, this would give a substantial boost to lattice 
QCD. Therefore, many groups have started to work on 
this problem. Here we present our first results using some 
interpolators which have not been tested on the lattice 
yet. 

In our calculations we perform a qualitative study using 
different types of spin- i operators with the quantum num- 
bers of the 0'^. We compute all cross correlators and use 
the variational method |33ll34j to extract the lowest lying 
eigenvalues. These are used to create effective mass plots 
for a comparison to the N-K scattering state which we 
computed separately on the same lattice. 



2 Details of the calculation 



We considered the following interpolating fields as basis 
for our correlation matrix: 



Currents suggested by Csikor/Fodor 



Iq/I = iabc[UaCj5db]{UeSeij5dc T ^ d)} (1) 



Currents suggested by Sasaki 



6>\- = e.abc<^aef(-bgh{ulCdf){ulC"fzdh)Cs'^ (2) 



f^abceaefebghiu'^ Cj5df){u^ Cj5j^dh)Cs'^ (3) 



Ot^f^ = iabceaefebghiu'^Cdf){u'^Cj5-ff_,dh)Cs'^ (4) 



Currents which are suggestions by L. Ya. Glozman j^Sl, 
however, using only s-wave quarks instead of a mixture 
of s-wave and p-wave quarks in jsj and (jSJ : 



Ifj. — {SaeSbg + SbeSag)igcd 
UaCub 

V2 



(5) 



X I :^{ulCdb + d^Cub) 



di Cdb 

X I j=^{u^{C^,)d, + dnCiM \CsJ 
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I = {^aehg + 5be5ag)^gcd (6) 

-^{ul{Cj^,j5)db + d^(C7^75)wb) 

Other currents related to suggestions by L. Ya. Gloz- 
man using only s-wave quarks and a factor {Sae^bg — 
SbeSag) iustcad the factor {SaeSbg + Sbe^ag) in (jZJ and 



Ifj. — {SaeSbg — SbeSag)egcd 

X [u^{Cj5)dd-d'^{CjMCs'^ 



Iv = {5ae5bg — SbeSag)£gcd 

X ( -7n{ui{Cof,y)db + d^{CcFf,y)ub) 



(7) 



(8) 



di iC'af,^)db 
dI{Cj^)dd 



X I ^iuJiCj^)dd + dJiCj,)ud) \Csl 



u^{Clpi)ud 



— While the interpolating fields Q to (0J were already 
used in other lattice calculations the interpolators ((SJl 
to (|HJl have not been tested on the lattice yet. 

In order to get only spin-^ pentaquarks, we have to 
project the spin. This is done using the spin projection op- 
erator for a Rarita-Schwinger field The corresponding 
spin-i state can be projected by applying the projection 
operator 



where 



(10) 

One can show that choosing the temporal components of 
the Lorentz indices at zero momentum is sufficient to get 
only spin-i contributions. 

Our baryon correlators are also projected to definite 
parity using the projection operator = ^(1 ± 74). We 
obtain two matrices of correlators: 

C^(t) = Z±e-^^^ + Zr^e-^-^^- , (11) 
when we project with and 

C^{t) = -Z^e-'^- -Z^e-^-^^\ (12) 



when we use P as shown in |87j . These two matrices are 
combined to 

C{t)^l{C+it)-C-{T-t)) , (13) 

to improve the statistics. This in combination with our 
quarks having antiperiodic boundary conditions in tem- 
poral direction allows us to separate the parity channels. 
The parity channel of the 63+ is not known. There are 
conflicting theoretical predictions and no conclusive ex- 
perimental data, so we have to look at both channels. 
For the interpolators lO to (O p-wave quarks^ would be 
required. Since we do not have p-wave sources we had to 
adjust the color structure in the interpolators {Tj) and (|SJ) 
to obtain a signal at all^. By doing so one immediately 
sees that the interpolator Q becomes, up to a constant 
factor, the same as the operator (PJ. Performing a small 
test simulation on a few configurations, we find that the 
interpolator ^ is numerically the same as interpolator 
(01. Therefore we exclude the interpolators (0) and Q 
from our analysis. 

The interpolator 10 contains two diquarks with / = and 
so it has isospin / = 0. In contrast the interpolators 
(01 and ((HJ are linear combinations of two diquarks with 
I — 1. Thus they are a mixture of isospin / = and 1 = 2 
states. 

We also have to exclude the interpolators (^ because the 
cross correlations between them and the other interpolat- 
ing fields are more complex and thus more demanding in 
computation. 

Then we use the remaining five interpolators to calculate 
a cross correlation matrix Cij(t) which is then inserted 
into the generalized eigenvalue problem 

a,{t)vf = x^''\t)a,{to)vf. (14) 

The solutions of this equation behave like 
(X exp 

These eigenvalues are used to compute effective masses 
according to 



"^e//(^) = In 



m 

A(t + 1) 



(16) 



Ordering the five eigenvalues according to their absolute 
value the largest eigenvalue in the positive parity channel 
should give the 0"*" mass if is a positive parity particle. 
The second largest eigenvalue in the negative parity chan- 
nel should give the mass, where the largest eigenvalue 
corresponds to the N-K scattering state at rest. 
In our quenched calculation we use the Chirally Improved 
Dirac operator |38II39| . It is an approximate solution of 
the Ginsparg- Wilson equation j4Q|, with good chiral be- 
havior |41II42II43) . The gauge fields are generated with the 

^ Quarks with an orbital p-wave excitation in their wave 
functions. 

^ We changed the color factor {SaeSbg + SteSag) into {SaeSbg — 

SbeSag)- 



Christian Hagen, Dieter Hierl, Andreas Schafer: Search for the ©^(1540) in lattice QCD 



3 



size X T 


12^ X 24 


P 


7.90 


a [fm] 


0.148 


a-^ [MeV] 


1333 


L [fm] 


f« 1.8 


#conf N 


100 




0.04, 0.05, 0.06, 0.08, 


quark masses aruq 


0.10, 0.12, 0.16, 0.20 


smearing parameters: 


n = 18, K = 0.210 


s-quark mass anris 


0.0888(17) 



Table 1. Parameters of our calculations. 



Liischer-Weisz gauge action |44II45| at /9 = 7.90. The cor- 
responding value of the lattice spacing is a = 0.148 fm as 
determined from the Sommer parameter in |46| . For the 
quark sources we use a Gaussian-type distribution gen- 
erated with Jacobi smearing |47i»48) . The error bars are 
computed using the jackknife method. The s-quark mass 
is determined from the Kaon mass. The parameters of our 
calculation are collected in Tabled 



3 Results 

The results of our calculations are shown in Fig. ^ where 
we plot the effective masses of the two lowest lying states 
of both parity channels obtained with the cross-correlation 
technique. These states are approaching a possible plateau 
very slowly as we expected, since states consisting of five 
quarks are very complicated and therefore should contain 
a large number of excited states which have to die out 
before the effective mass reaches a plateau. We use in ad- 
dition to the cross-correlation technique Jacobi smeared 
Gaussian quark sources for all our quarks to improve the 
signal for the lowest lying states. 

The lower dashed horizontal line in the negative parity 
channel is the sum of the nucleon and kaon mass at rest 
in the ground state obtained from a separate calculation 
on the same lattice. Since we project the final state to zero 
momentum a scattering state can also be a two particle 
state where the two particles have the same but antipar- 
allel momentum, i.e. pjv — —pK- We use the relativistic 
£'-p-relation to calculate the energy of such states, 

E = ^Jp^+m% + ^p'^+m]^, (17) 

where the smallest momentum is 2tt/L « 700 MeV on our 
lattice. In Fig. ^ this energy is represented by the upper 
horizontal line. 

In contrast to similar calculations we performed for or- 
dinary mesons and baryons we do not find clean mass 
plateaus in this case. This implies also that we cannot 
perform any systematic chiral extrapolation. The signals 
we get are compatible with N-K continuum states, but 
do not unambiguously identify them. For the latter we 
would have needed results for different volumes, which 
would have been very expensive. Furthermore, we are not 



interested in identifying such a continuum state, but only 
want to know whether there is any indication for a novel 
bound state. 

In the negative parity channel we find effective mass pla- 
teaus which are consistent with the lowest N-K scattering 
states. We find that the second state is noisy but within 
errors consistent with the energy in Eq. 1)17(1 . Therefore, 
it is most likely that we do not observe a 6*+ state in the 
negative parity channel. Naturally, if the were broad, 
which implied that it would mix strongly with the contin- 
uum our conclusion would be weakened. However, for the 
Q'^ this possibility is excluded experimentally. 
From the composition of the eigenvectors, shown in Fig. 
|5J we conclude that the two lowest lying states are really 
independent of each other. One can also see that some of 
the interpolators which have not been used in previous 
studies actually give large contributions to the low lying 
states. 

In the positive parity channel one expects to find either a 
bound 0+ or an excited N-K scattering state. For such 
an excited state there are several possibilities, e.g., N*-K, 
or N-K with relative angular momentum, and so on. 
On the positive parity side, we also show the two lowest 
lying states obtained from our calculation. Both of them 
are too heavy to have anything to do with the 0"*". They 
probably correspond to excited N-K scattering states. If 
there were a signal belonging to the it is supposed to 
lie below the solid line assuming that the chiral extrapo- 
lation of the 0'^ does not lead to dramatic effects below 
our smallest quark mass. (A serious chiral extrapolation 
would require cleaner mass plateaus as noted above.) 

4 Conclusion 

In this article we present the results of a pilot study for 
searches on the lattice using different types of oper- 
ators. We find that for negative parity our results are in 
good agreement with a N-K scattering state in the ground 
state and a quite noisy signal for the first excited state. 
For positive parity we find states which are typically more 
than 500 MeV heavier than the 0~^ and thus not compat- 
ible with a 6>+ mass of 1540 MeV. ^ 

In our calculation we only used one lattice volume. If we 
had found a 6)+ candidate, for the final distinction be- 
tween a bound or a scattering state we would have had 
to study the volume dependence. As it stands, there is no 
motivation to do so. 

Up to now our calculation do not show any hints for a 0^ 
in the quenched approximation with chiral fermions using 
in addition to standard ones the new interpolators Q-llHl)- 
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Fig. 1. Results from cross-correlation of our final set of interpolating fields, i.e., JIJ, 0, 0, Q and Q. We show effective 
masses for the two lowest- lying eigenvalues computed according to Eq. 1161 . The dashed line represents Mjv -f Mk obtained 
from a separate calculation on the same lattice. The solid line is the energy for the smallest momentum calculated according to 
Eq. 1171 . For all our quarks we use Jacobi smeared Gaussian sources. 
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Fig. 2. Plot of the eigenvector components for amq — 0.08 
which corresponds to the pion mass m,r = 0.66GeV. The left 
hand side plot shows the eigenvector components of the neg- 
ative parity ground state, the right hand side plot is for the 
negative parity 1st excited state. The eigenvectors are obvi- 
ously stable and different. Thus, we conclude that we observe 
to independent states. 
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